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Recent experiments on Bose-Einstein condensates in optical cavities have observed a coherent 
state of the matter-light system — superradiance. The nature of these experiments demands con- 
sideration of collective dynamics. Including cavity leakage and the back-reaction of the cavity field 
on the condensate, we find a rich phase diagram including a variety of multi-phase co-existence 
regions, and a regime of persistent optomechanical oscillations. Proximity to some of the phase 
boundaries results in critical slowing down of the decay of many-body oscillations, which can be 
enhanced by large cavity loss. 

PACS numbers: 37.30.-K, 42.50.Pq 



The tremendous advances in preparing Bose-Einstein 
condensates (BEC) in optical cavities have opened new 
frontiers combining cold atoms and quantum optics. It 
is now possible to enter the strongly coupled regime of 
cavity quantum electrodynamics (QED) [l|, [i[ in which 
atoms exchange photons many times before spontaneous 
emission and cavity losses set in. The inherent cav- 
ity leakage also provides a valuable window on these 
quantum many-body systems. In particular, it allows 
for in situ non-demolition measurements of condensate 
properties via optical transmission Q). The strong 
matter-light coupling also supports collective dynamics 
and back-reaction effects, stimulating new directions in 
cavity optomechanics [H, @] and self-organised atomic en- 
sembles 0413. 

More recently, these capabilities have been elevated 
through observation of the superradiance transition in 
BECs [13, 15 1 . The atom mediated coupling between 



a transverse pump field and a cavity mode leads to a 
realisation of the Dicke model 16Hl9f. in which atomic 
momenta play the role of spin states; see Fig. [TJ A sig- 
nificant merit of this approach is that the energy split- 
ting of the two-level systems is small enough so that the 
Dickc superradiance transition may be realised with light 
at optical frequencies Il5j |. These experiments are a 
landmark in the study of quantum phase transitions in- 
volving spins, and offer exciting and unique prospects 
for exploring their static and dynamic properties. In- 
deed, the time-dependent nature of these measurements 
demands consideration of collective dynamics. 

Motivated by these developments we investigate the 
collective dynamics of BECs in optical cavities. Our 
two primary goals are to establish the generic behaviour, 
and to focus on the precise experimental realisation in 
Ref. [HI . We obtain a surprisingly rich phase diagram for 
a broad range of parameters, and find distinct regimes of 
dynamical behaviour, including several regions of multi- 
phase co-existence, and regions of persistent optomechan- 
ical oscillations. For recent theoretical work see Ref. 1201. 
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FIG. 1. Experimental setup [15||- (a) BEC in a transversely 
pumped cavity with pumping frequency, u p , and strength fl p , 
single-atom cavity coupling go, atomic transition frequency 
u) a , cavity frequency ui c , and cavity decay rate, n. (b) Energy 
levels and pumping scheme showing the two-level splitting, 
luo = 2uj r , in the effective Dicke model, where uj r = k 2 /2m is 



the recoil energy. 



with N ~ 10 5 atoms in an optical cavity with a transverse 
pumping laser; see Fig. [1] The excited atoms may re- 
emit photons either along or transverse to the cavity axis. 
This process couples the zero momentum atomic ground 
state, \p x ,Pz) — |0, 0), to the symmetric superpositions 
| ± k,±k), with an additional photon momentum along 
the cavity or pump directions. This yields an effective 
two-level system or "spin", where the splitting, ujq, is 
twice the atomic recoil energy, ui r = k 2 /2m. One obtains 
an effective Dicke model for collective spins, S, of length 
N/2, coupled to radiation ip 14, 15| 



The experiments in Ref. 15| consist of a 87 Rb BEC 



H = unpty + uj S z + US z ip^ip + g(^S~ + ipS + ) 

+ g'(^S+ + if>S-), (1) 

where, u) — u) c — uip + NUq(1 + A4)/2, ojq = 2u) r , U = 
UoM, M is a matrix element of order unity, and Uq = 
go/(up — Lu a ) encodes the back-reaction of the cavity light 
field on the BEC. The model includes both co-rotating 
and counter-rotating matter-light couplings, denoted g 
and g' . In the experiment g = g' = goQ p /(co p — Lo a ) [HI]. 

To describe the dynamics of the matter-light system 
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(TTJ we construct the Heisenbcrg equations of motion 

S- = -i(uo + U^i>)S~ + 2i(gi/j + g'^)S z , 
S z = -igi}S+ + ig^S' + ig'^S' - ig'^S + , (2) 
ip = - [as + i(u + US Z )} V> - igS- - ig'S + , 

where S = S x US,, , k is the cavity loss rate, and we 
neglect atom loss [Iff. Various limits of these equations 
have been explored in different contexts. For k = g' = 
they describe fermionic pairing, where _0 is the Feshbach 
resonant closed state molecular field [21j. This regime 
also arises for polariton condensates and phase-locking 
of oscillators [22J. More recently, for g = g' , they have 
emerged in an elegant proposal for realising the Dicke 
model 14 1 . As we will see, solutions of the more general 
equations strongly influence g = g' dynamics. 

In order to anchor the complete phase diagram, we 
start with {7 = and consider {7^0 below. Numerical 
solution of equations and the arguments below, yield 
the rich phase diagram in Fig. [2j where the phases indi- 
cate stable attr actors of the long time dynamics. Four 
distinct phases exist corresponding to all spins down and 
no photons (JJ-), all spins up and no photons (ff), a super- 
radiant phase with photons (SR), and co-existence of the 
supperradiant and down attractors; see S1-S4 in Fig. [2] 
Such co-existence, or bistability, is related to the ob- 
served optomechanical oscillations in a different system, 
where matter-light coupling is only through the back- 
reaction 0, [23[. In spite of the cavity decay rate, as, 
which may be large, the counter-rotating terms stabilise 
superradiant steady states. Indeed, dropping the deriva- 
tives in equation ([2]) renders algebraic equations, and the 
determinantal condition for non-trivial solutions (tp ^ 0) 
yields 



S z 



-0JLU (g 2 +g' 2 )± v/(2u;w <?<?') 2 
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2(g 2 - g' 2 ) 2 



(3) 

The conditions for real physical solutions yield the blue 
phase boundaries shown in Fig. [5](a). Setting S z = —N/2 
in equation ([3]) yields the "upper" boundary shown in 
Fig. [2](a) [24j. The vanishing of the square root yields 
the "lower" boundary, g' = g^f ct-/a +1 where a± = 
\/w 2 + k 2 ± cj, delineating the onset of co-existence. In 
order to identify the green phase boundary in Fig.[5](a) it 
is necessary to consider the stability of the steady states. 
We consider fluctuations about an arbitrary configura- 
tion, S = So + <5S, V = "00 + ty, with frequency v. 
Instability occurs if Im(V) > 0, and this yields the criti- 
cal line g 1 = g^J p+/ /?- shown in Fig. [5](a), where (3± = 
(uj ± uq) 2 + k 2 . This separates the stable normal state JJ- 
from the stable inverted state ff. For the chosen param- 
eters this gives g'/g = 1.0043, very close to unity. The 
dynamics at g = g' may thus be strongly influenced by 
proximity to this phase boundary. The parameters used 
in Fig.[5]follow the hierarchy w,k> ffV^V" 3> wo 5 m which 




0.0 0.5 1.0 1.5 

gTN (MHz) 



FIG. 2. Dynamical phase diagram for (7 = and correspond- 
ing spin trajectories on the Bloch sphere, (a) Dynamical 
phase diagram for parameters to = 20MHz, too = 0.05MHz, 
ft = 8.1MHz taken from Ref. [l5(|. showing the stable at- 
tractors of the nonlinear dynamics for U = 0. The phases 
are: ij. all spins down S z = —N/2 and no photons, -ft all 
spins up S z = N/2 and no photons, a non-trivial magne- 
tised superradiant state with %j> ^ 0, and a co-existence re- 
gion emanating from a tricritical point •. The separatrix 
g'/g = \//3+//3- = 1.0043 is close to but distinct from 
unity, (b) Small cavity loss regime with k = lKHz show- 
ing the evolution towards the superradiance transition at 
g + g' — sJZJuJa/N in the equilibrium Dicke model [l^ - H^ ]. 
The Bloch spheres S1-S4 show the stable •, unstable o, and 
hyperbolic x fixed points (steady states) as well as character- 
istic trajectories in each of the phases. Examples of the time 
evolution for g = g' are given in Fig. [3] 



the photon decay rate, k = 8.1MHz, is much greater 
than the level spacing, ujq = 0.046MHz [l5j]. In this 
limit one obtains a characteristic decay rate for the col- 
lective many-body oscillations, Im(j/) = —ku>q/(k 2 + uj 2 ), 
as verified in Fig. [3] (b). Notably, in the limit k — > oo, 
corresponding to a large cavity loss rate, this results 
in Im(/y) — > 0, or slow decay of the collective oscilla- 
tions. This may be understood as critical slowing down 
[25I ]. Further insight into this k — > 00 dynamics may 
be gained by adiabatic elimination of the fast photon 
field, i> = -[i(g + g')S x + (g - g')S y }/{K + iuj), to de- 
rive an effective equation of motion for the classical spins 
S = {S,i7}-rSx(Sxz). Hcrcff = ujqS.-A+S^-AS 2 
is the Lipkin-Meshkov-Glick Hamiltonian [26|, |27| , with 
A± = 



■{g±g'f andT 



r( S J -g 2 )- The ad- 



K 2 +UJ 2 V y & > K 2 -\-U) 2 

ditional term takes the form of damping in the Landau- 
Lifshitz-Gilbcrt equations [28|. Depending on the sign 
of r this favours spin alignment either parallel or anti- 
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FIG. 3. Time evolution of the effective spin state and cav- 
ity photon population, (a) Evolution of S z and (b) photon 
number for g-/N = g 'y/N = 0.791MHz and [7 = 0. The 
long time behaviour shows the exponential envelope \ip\ 2 = 
\ipo\ 2 ± Ae lm{v)t (dashed lines) where A is a non- universal 
amplitude dependant on the initial conditions, and the decay 
rate \m(v) — — kujq/(k 2 + oj 2 ). (c) For ujo <C to, the long time 
oscillation frequency is well described by the perturbative re- 
sult Re{v) = u \S\/S z + 5, where 8 = Aug 2 S 2 z /\S\{K 2 + ui 2 ) 
is a small correction to leading term of order ujq. The short 
and intermediate time dynamics can be strongly affected by 
the existence of additional stable or unstable fixed points. 



parallel to the z-axis. The sign change at g = g' is 
consistent with the n — > oo limit of the phase bound- 
ary, g' = g^J (3 + / ft-, which separates the JJ- and ft steady 
states. It is interesting to contrast the emergence of intc- 
grable dynamics for g = g' and k oo, with the chaotic 
behaviour when g = g' and k = [l9(. Moreover, for 
g =/= g' the dynamics is non-Hamiltonian. 

Having discussed the dynamics of the model ([1]) for 
U = 0, let us now consider U ^ 0. In order to make close 
contact with the experimental realisation in Ref. (ljj we 
henceforth set g = g'. In Fig. |U(a) we present the dy- 
namical phase diagram as a function of U. The entire 
topology may gleaned analytically from the steady state 
solutions of equation @. These reveal two classes of su- 
perradiant solutions incorporating both U and k. The 
first class has a non- vanishing photon population 
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whe 



re ui = oj 



US Z , and 




g 2 (4uj 2 -U 2 N 2 )-Uuj k 2 



U 2 (uj U + 4g 2 



S y = 0, (5) 



and S x is determined by the normalisation of the spins. 
Physical solutions require < N/2. In the limit 

k, U — > we recover the results of equilibrium super- 
radiance IB" 13; an d f° r U — they reduce to those 
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FIG. 4. Dynamical phase diagram as a function of U and long 
time steady states, (a) Dynamical phase diagram of model 
(P) with g = g' and parameters w = 20MHz, lj = 0.05MHz, 
K = 8.fMHz taken from Ref. [lfj. The blue, red and gold 
critical boundaries correspond to equations Q. (b) Magnified 
portion showing the appearance of a bistable superradiant 
phase (2SR) corresponding to co-existence of both roots of 
equation ([5]). (c,d) Trajectory along the dashed-dotted line 
showing the comparison between steady state solutions and 
numerical integration of the equations of motion at 360ms. 
The region to the right of the blue asymptote corresponds 
to a limit cycle. For each value of U we take a variety of 
initial conditions with -0 = 1 and S uniformly distributed 
over the Bloch sphere, (e) Magnified portion of the bistable 
superradiant phase (2SR) showing the agreement between the 
steady states and numerical integration. 



of Ref. [141 ]. For sufficiently large negative U it is pos- 
sible for equation (|5|) to develop unphysical complex 
roots. In this case one may satisfy equations ([2]) with 
u} =cj + U\4'\ 2 = 0, u = uj + US Z = 0, and 







U 



S z 



-v- < 6 » 



where S y is determined by normalisation. Physical so- 
lutions have S 2 + S 2 < N 2 /4. In general these distinct 
solutions are connected for g =/= (/, so we do not distin- 
guish them in Fig. H](a). Nonetheless, it is important to 
keep track of them for analytical work when g = g' ■ Fig- 
ure S](a) consists of three phase boundaries correspond- 
ing to instability of 4 (blue), instability of ft (red), and 
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FIG. 5. Dynamical phase diagram as a function of u for the 
experimental parameters used in Ref. [l5|j. Dynamical phase 
diagram in the gy/N versus ui plane for UN — — 40MHz, 
where the blue, red and gold phase boundaries are given by 
equation {?} and correspond to those in Fig. |4](a). The thick 
blue line is the boundary of stability of the JJ. state that would 
be seen on increasing g as in Ref. [IjJ. Inset: Phase diagram 
re-plotted as a function of g 2 N for comparison with Fig. 5 of 
Ref. [II. 



existence of the second-type superradiant phase (gold): 



value of U. We see that the supcrradiance boundary is ac- 
companied by several regions of multi-phase co-existence. 
It would extremely interesting to investigate this exper- 
imentally. The inset shows the same data shifted and 
rescaled for comparison with Fig. 5 of Ref. [l5| . 

In summary, we have discussed the collective dynam- 
ics of BECs in optical cavities. We obtain a rich phase 
diagram with different regimes of dynamical behaviour, 
including several regions of multi-phase co-existence and 
the slow decay of many-body oscillations. Amongst our 
findings is a regime of persistent optomechanical oscil- 
lations described by a damped driven pendulum. Given 
the strong interest in cavity optomechanics 0, Q this may 
be a profitable region to explore experimentally. Further 
directions include the impact of cavity axis pumping [29| 
and photon correlations. Experiments in which the cou- 
pling g is quenched through the phase boundaries may 
help explore this rich dynamics. 

We are extremely grateful to K. Baumann, F. Bren- 
necke, T. Esslinger and M. Kohl for illuminating dis- 
cussions. MJB and JK acknowledge ETH Zurich, G. 
Blatter, S. Schmidt, and H. Tiireci for hospitality and 
interactions. MJB and BDS acknowledge EPSRC grant 
no. EP/E018130/1. JK acknowledges EPSRC grant no. 
EP/G004714/1. 
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where lojj = UN/2. Instability of the normal state g^ 
has also been considered for thermal clouds in a ring 
cavity The result for g* delimits the region, both 
for equation ([5]) and equation (|6]), to have real, phys- 
ical solutions. All three of these boundaries intersect 
at U 
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; 2 , g = \/—(jJoU /4, as shown in 
Fig. |4](a). Upon increasing g one finds a phase where 
two distinct first- type superradiant solutions co-exist; see 
Fig. |4](b). This is borne out in Fig. H](c,d,e), where we 
compare the steady states with direct integration of equa- 
tions ([2]) along the dashed-dottcd line in Fig. |4](a). We 
integrate over a period of 360 ms to eliminate the tran- 
sitory effects of critical slowing down discussed earlier. 
In addition, Fig. [4] contains several regions involving co- 
existence of superradiant and non-superradiant phases. 

The scattered points in Fig.@](d), beyond the J| bound- 
ary at U — 2uj/N, correspond to limit cycles rather than 
steady states. Here S z = —u/U and -0 is purely imagi- 
nary. Writing S~ = re~ l8 , where r = -J ' N 2 /4 — uj 2 /U 2 , 
yields d t 9 = ujq + U\ijj\ 2 , and (d t + n)ip = —2igrcos(9), 
with limit cycle behaviour. For n ^> luq + U\ip\ 2 , these 
describe a damped driven pendulum. 

Having confirmed the overall phase diagram in 
Fig. |4](a) as a function of U, let us finally focus on the 
specific value UN = -40MHz used in Ref. [3]. In Fig. [5] 
we plot the phase diagram as a function of uj for this fixed 
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